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Abstract. We present in this paper a set of routines con- 
structed to compute the rank of a matrix Lie algebra and also 
to determine a Cartan subalgebra from a given list of ele- 
ments. 

■ Introduction 

It is clear the usefullness of the Cartan subalgebras in the Lie algebra 
context (see the references). We have worked in a previous paper in the 
links between the Lie algebras and particle physics. In that work we needed 
a routine to construct Cartan subalgebras (in order to define observables) 
and, as a consequence, the rank of a Lie algebra. In this paper we present all 
the routines defined for that purpose. 

We will show some examples with clasical Lie algebras (like so(n), sl(n)) 
and with two models in particle physics (split g2 y so(4, 4)®so(2, 2)). The 
theoretic result on which the rountine is based is the following: 

■ Teorema. 

Let £ be a real semisimple matrix Lie algebra finite -dimensional Let H be 
a subalgebra of L. Then H is a Cartan subalgebra if and only if it is 
abelian maximal and for all x and y in H, [x, y' ] = . 

The routine proceeds as follows: first of all a list of elements of the algebra 
is introduced. Once the necessary conditions are verified, its maximality is 
analyzed. If it is maximal, they form a Cartan subalgebra and its cardinal is 
the rank of the Lie algebra, but if they are not, the routine solve the system 
for the possible elements that can extend the set in order to form a Cartan 
subalgebra. This process is repeated until the subalgebra is maximal. 

■ Auxiliary commands 

The following command computes the intersection of two subspaces. : 



Rank of a matrix Lie algebra 



3 



«5i_,j_ := If [i == j, 1, 0] ; 
Cap[m_, n_] : = Module [{} , 

Length[m] Length[n] 

genl := ^ Xi m|[i]| ; gen2 := ^ y±n|[i]|; 

i=l i=l 

in = genl // . ToRules [Reduce [genl == gen2] ] ; 

varia : = Variables [in] ; num : = Length [varia] ; 

rules := Table [varialjj -» <5i, j , {i, 1, num}, {j, 1, num}]; 

Table [in//. rules[[i]], {i, 1, num}]]; 

and of course the Lie bracket: 

c[x_, y_] : = x.y - y.x; 

We start from a matrix Lie algebra £ with basis B constituted by the matri- 
ces {t>j} n j _ l (previously defined). The future basis of the Cartan subalgebra 
"H will be denoted by (v,} r , where r is the rank. The next function, 
KerAd[m,b,n] where m is a matrix, b is the notation for the basis and n is 
the dimension of the algebra, computes de kernel of the linear map: 

ad m : X £ 
x i — > [m, x] 

where [m, x] is the Lie bracket in £, that is to say, [m,x]=mx-xm: 
KerAd [m_, b_, n_] : = Module [ { } , 

n 

X = ^a ± bi ; alia = Table [a t , {t, 1, n} ] ; 

i=l 

Matr = { } ; 

Do[AppendTo[Matr, alia//. Solve [X == c [m, bi ] ] PI ] , {1, 
Matr = Transpose [Matr] ; NullSpace [Matr] ] 

The next module defines Indepen[v,r], which checks if the elements are 
linearly independent: 

Indepen[v_, r_] : = Module [{fal = 1, indep = { } } , 

r 

indep = ToRules [Reduce [^d k v k == e] ] ; 

k=l 

Do [If [ ( (d k // . indep) = = = 0) == False, fal = 0] , {k, 1, : 

The next one, Co[v,r], check if the set is abelian: 



Pablo Alberca -Cdndido Martin 



Co [ v_, r_] : = Module [ { f al = 1 } , 

Do [If [c[v ± , Vj ] * e, fal = 0] , {i, 1, r}, { j, i + 1, r}] ; fa 
The module Star[v,r] check if the set {v,} r , satisfies the condition 

7-1 

[x, x'] = for all x: 

Star[v_, r_] : = Module [ {fal = 1} , 
Do[If[c[vi, Transpose [vj ] ] # ©, fal = 0] , 
{i, 1, r}, {j, i, r}]; fal] 

Now, the command is Maximal[v,r,b,n] which establishes when an abelian 
set is maximal. : 



Maximal [v_, r_, b_, n_] : = Module [ {fal = 1} , 

ec = {c[X, Transpose [X] ] ==©}; 

Do[ 

AppendTo[ec, c[v k , X] ==6], {k, 1, r}]; 
Do [AppendTo [ec, 
c [Transpose [v k ] , X] ==6], {k, 1, r}]; 

Do [AppendTo [ec, 
Scala[X, v k ] == 0] , {k, 1, r}] ; 
rulesol = Reduce [ec] ; 

If [Length [rulesol] # n, fal = 0] ; 
fal] 

At last, we have the module Addone[r,v,b,n] which adds, inside the con- 
struction process of a Cartan subalgebra, a new element in each step: 

Addone [r_, v_, b_, n_] : = Module [ { } , 
verif :=c[Y, Transpose [Y] ] == 9; 
Do [verif = verif && c [Y, Vj ] ==6 && 

c[Y, Transpose [vj ] ] == 6 &&Scala[Y 

{j, r-1}]; 
ta = Y // . {ToRules [Reduce [verif ] ] } PJ ; 

var = Variables [ta] ; 
num = Length [var] ; 

rul = Table [var [[jj -» 5i, j , {i, num}, {j, num}]; 

kk = 1 ; 

While [ (ta // . rul[[kk]]) == 6 , kk = kk + 1] ; 
v r = ta //. rullkk]]] 



■ Main routine 

It is time to present the main routine. We need first to introduce a basis of 
the particular Lie algebra we are interested in. We are going to present a 
comand for each algebra but it is important to say that, of course, the routine 
remains the same. 



Rank of a matrix Lie algebra 



5 



■ Lie algebra so(n) 

We compute first a basis of the Lie algebra: 

BSo[t_] : = Module [{kk = 1} , 
1 

n = — (t - 1) t; = Table[0, {i, 1, t}, { j, 1, t}] ; 
2 

Do[m = e; m[[p, qfl = 1; F p , q = m, {q, 1, t}, {p, 1, t}] ; 
Do[J p , q = F p , q - F q , p , {p, 1, t}, {q, 1, t}] ; kk = 1; 
Do[Do[b kk = J p , q ; kk = kk + 1, {q, p + 1, t}] , {p, 1, t - 1}] ; 

t t n 

Scala[x_, y_] :=^^x[[i, jj y[[i, jj ; X = £ a± b± ; ] 

i=l j=l i=l 

and the the main routine is: 

RankOf So [t_, dat_] : = Module [ { } , 

r = Length [dat] ; 

Do[v k =dat[[k]], {k, 1, r}]; 

If [lndepen[v, r] ==1&&Co[v, r] ==l&SStar[v, r] = = 
pset = KerAd[vi , b, n] ; If [r > 1, 

Do [pset = Cap [pset, KerAd[vi, b, n] 

Length [pset] ^ 

Y= ^ f k ^psetlk, i]] b ± ; 

While [Maximal [v, r, b, n] ==0, r = r + 1; Addone[r, v, b, i 
Print ["The rank is ", r, " and a Cartan subalgebr 

Do [comb = Solve [v q == X] ; 

n 

Print["v" q , " = ", ^ (a p //. combll]] ) * b p ] , {q, 1, i 
Print [False] ] ] 

The fundamental line is 

While [Maximal [v, r, b, n] == 0, 
r = r + 1; Addone [r, v, b, n] ] 

Let's look an example: 
BSo[6] 
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and then a generic element is 



X 



( 


ai 


02 


fl3 


a 4 


a 5 s 


-a.\ 





a 6 


fl 7 


a 8 


a 9 


-a 2 


-a 6 





aio 


a n 


ai2 


-fl3 


—a-j 


-aio 





an 


«14 


-04 


-as 


-a n 


-an 





ai5 


> -fl5 


—at) 


-an 


-an 


-1215 


, 



We make then 



RankOfSo[6, {bi}] 

The rank is 3, and a Cartan subalgebra: 



v, =b, 
v 3 = b l0 

where we have also shown the time it took. We have then that the rank is 3 
and a Cartan subalgebra is constituted by the matrices 



{ {bi , bio i bis } } 

II 1 0] 

-10 







,1 0, 



(0 0] 



1 
0-1000 



,0 0, 



(° 














°1 







































































1 


,0 











-1 


o, 



We can start form two elements, b\ y b l0 , and we have 
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RankOfSo[6, {bi , bio } ] 

The rank is 3, and a Cartan subalgebra: 

vi =bi 

v 2 = b m 

v 3 =b l5 
Another example could be 

BSo[9] 

and the a generic element is 



X 








a 2 


a 3 


fl4 


a 5 


a 6 


a 7 


at s 


fl] 





ag 


aw 


a n 


an 


an 


an 


a\s 


a 2 


-Og 





«16 


a [7 


an 


«19 


a 20 


«21 


a 3 


-«io 


_ «16 





a 22 


a 2 3 


a 2 4 


a 2 s 


«26 


fl4 


-a n 


-an 


-a 22 





a 21 


«28 


a 2 g 


«30 


as 


-an 


-an, 


-an 


-an 





«31 


a 32 


«33 


a 6 


-an 


-a\g 


-a 2 4 


-a 2 z 


-a 3l 





«34 


fl 35 


a-] 


~ a u 


~ a 20 


~a 2 s 


~a 2 g 


~a 32 


-fl 34 





«36 


at 


-a\s 


-a 2 i 


-a 2 e 


_ «30 


-a 33 


-a 35 


-a 3 6 


o , 



We can make then 

Rank0fSo[9, {b x }] 

The rank is 4, and a Cartan subalgebra: 

Vi =b x 

V 2 = ^36 
V 3 = b 3l 

v 4 = b 22 
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As a final example, we can use the Lie algebra so(32), recently used in 



32 

string theory (dimension=^- 31 = 496): 



BSo[32] 

and a generic element 



""93 ""120 ~"\ 

""94 -"121 

""95 ""122 

""96 -"123 

""97 ""124 



-"m 

""129 
""130 
""131 
""132 



""197 
""198 
""199 
""200 



-"203 
-"204 
""205 
""206 
""207 



-"221 
-"222 



'"226 
'"227 
""228 
""229 
""230 



"222 
"244 



"200 
"223 
"245 
"266 



"126 
"152 
"177 
"201 

"224 
"246 
"267 



-"269 
""270 
""271 
""272 
""273 



"178 
"202 
"225 
"247 
"268 
"288 
"307 



""2S9 

-"290 

""291 
""292 
""293 



""309 
""310 

-"311 
-"312 



"179 

"203 
"226 
"248 
"269 
"2H9 



""344 
""345 
""346 
""347 



"204 
"227 
"249 
"270 
"290 



""361 

""362 

-"363 



"250 
"271 
"291 
"310 

°328 

"345 
"361 


""377 
""378 



"182 
"206 
"229 
"251 
"272 



"329 

"346 
"362 
"377 



"207 
"230 
"252 
"273 



"347 
"363 

"378 
"392 



-"138 
""139 
-"140 



""211 
""212 
""213 
""214 
""215 
""216 



'"234 
""235 
-"236 
-"237 
-"238 
-"239 



""257 

-"258 

""259 
-"260 
-"261 



""299 
""300 
""301 
""302 



""3 IX 
""319 
""320 



""336 ""353 
""337 ""354 
""338 ""355 



" a 3S2 
""383 

""3Si4 



""396 ""409 

""397 ""410 

""398 ""411 

""399 ""412 

-0400 ""413 



■"401 



"414 



""144 
""145 



""219 
""220 



""264 
""265 



""305 
""306 



""324 
""325 



""359 
""360 



""375 
""376 



""390 
""391 



-"404 
""405 



'"417 
-"418 
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RankOfSo[32, {bi}] 

The rank is 16, and a Cartan subalgebra: 

vi =bi 

v 2 = b 496 

v 3 = b 49l 

v 4 = b 4S2 

v 5 = b 469 

v 6 = ^452 

v 7 = fc 431 

v 8 = b 406 

v 9 = b 311 

v 10 = &344 

v ll = ^307 

v 12 = ^266 

v 13 = ^221 

v 14 = ^172 

V [5 = b U 9 

Vl6 = b 62 



Lie algebra sl(n) 

The second example could be the Lie algebra sl(n). First of all we compute 
a basis as follows 
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Bsl[t_] : = Module [{k = 1} , n = t 2 - 1; 
6 = Table[0, {i, 1, t}, {j, 1, t}]; 

Do[m = 6; m[[p, qj = 1; F p , q = m, {q, 1, t}, {p, 1, t}] ; 
Do[Do[If [p * q, b k =F p , q ; k = k + 1] , {q, 1, t}], {p, 1, t}] 
Do[b k = Fi,i - F p , p ; k = k + 1, {p, 2, t}] ; 

t t n 

Scala[x_, y_] :=^^x[[i, jj y[[i, jj ; X = a± b ± ; ] 

i=l j=l i=l 

Then we have 

RankOfsl[t_, dat_] : = Module[{}, 

r = Length [dat] ; 

Do[v k =dat[[k]], {k, 1, r}]; 

If [lndepen[v, r] ==l&SCo[v, r] ==l&SStar[v, r] = = 
pset = KerAd[vi , b, n] ; If [r > 1, 

Do [pset = Cap [pset, KerAd[vi, b, n] 

Length [pset] ^ 

Y = V fk ^pset[[k, i]] b ± ; 

k=l 1 = 1 

While [Maximal [v, r, b, n] ==0, r = r + 1; Addone[r, v, b, i 
Print ["The rank is ", r, " and a Cartan subalgebr 

Do [comb = Solve [v q == X] ; 

n 

Print["v" q , " = ", ^ ( a P //• combll]]) *b p ], {q, 1, i 

Print [False] ] J 

We can have a basis of sl(3) by making 

Bsl[3] ; 

and then a generic element is 

X 



' a-/ + a% 


ai 


fl2 ' 


a 3 






> a 5 


a 6 


-«8 - 



If we take the first element of the basis, b\ , we obtain 
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RankOfsl[3, {b x }] 

False 

This is so becasuse 

c [bi , Transpose [bi ] ] 

1 1 Q\ 
0-10 

,0 0, 

and it must verify [x, x' ] = . Nevertheless, if we take b% : 

RankOfsl[3, {b 7 }] 

The rank is 2, and a Cartan subalgebra: 

v, =b 7 
v 2 =b s 

we get that the rank is 2 and a Cartan subalgebra is 
{ {b 7 , b 8 } } 



/ 


f 1 









( 1 





0^ 









-1 


















V 


,0 





o, 




,0 





-I) 


J 



In sl(5) we make first 

Bsl[5] ; 

with generic element 
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+ fl23 + «24 


a\ 


a 2 


fl 3 


a 4 N 


a 5 


-an 


a 6 


a 7 


a 8 


ag 


a w 


~a 2 2 


«ll 


«12 


«13 


a u 


«15 


-«23 


a 16 


a l7 


an 


«19 


<J20 


-#24 y 



Then we can make 



RankOfsl[5, {b 2 i , b 2 2 } ] 

The rank is 4, and a Cartan subalgebra: 



Vl =»21 
V 2 = 2 2 
v 3 = °23 ~~ ^24 

2 fo 2 l 2 fo 2 2 



V 4 



+ b n + b 7 . 



3 3 

to include these two elements in a Cartan subalgebra. 



Lie algebra g2 

We present now the Lie algebra g2 ■ A basis of this Lie algebra is 
Bg2 : = Module [ { } , 

n = 14; t = 8; e = Table[0, {i, 1, t}, { j, 1, t}] ; 
Do[m = 0; m[[p, qj = 1; F p , q = m, {q, 1, t} , {p, 1, t}] ; 
nnT.T = F - F /n 1 . 1-^ Sa 1 1-^ 1 • 



DO [— . — ~, v -xu — — r -p,q — "* r \.~ir —r — j / t, 

Do[J p , q =F p , q -F q , p , {p, 1, t}, {q, 1, t}]; 

bl = J2,3 + ^6,7 ; b 2 = J2,4 + J6,8 b3 = J 2; s + J7 ;4 ; 
b4 = J2,6 + ^8,4 ; bs = J2,7 + J4,5 ', ^>6 = Jz,8 + J4,6 ; 
b7 = J3,4 + ^7,8 ; bs = J3,S + J4,6 ^>9 = J3 ; 6 + J5,4 ; 
bin = J3,7 + J8,4 i bn = J3,8 + ^4,7 r bl2 = Js,6 + J7,8 



bl3 = Js,7 + J8,6 ; bl4 = J5,8 + ^6,7 ; 

t t n 

X!Zl x|[i ' 3l yli, jl; x = ^]a ± b i; ] 

i=lj=l 



Scala[x_, y_] 
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and then a generic element is, after making, 
Bg2 

as 

X 























"\ 








a [ 


a 2 




a 4 




a 6 





-a\ 





a 7 




ac, 


«10 


an 





-a 2 


-a 7 





fl 5 -fl 9 


a 6 + ag 


an -a 3 


-as, - aio 





-a 3 


-a & 


a 9 -a 5 





an 


an 


an 





—as, 


-at) 


-a 6 - a 8 


-an 





ai + 014 


a 2 - ai3 





-a 5 


-a 10 


a 3 — 


-a l3 


—a\ — a\s. 





a 7 +a 12 


.0 


-a 6 


-an 


O4 + O10 


-a l4 


fl!3 -a 2 


-a 7 - an 


, 



We define then 

RankOf g2 [dat_] : = Module [ { } , 

r = Length [dat] ; 

Do[v k =dat[[k]], {k, 1, r}]; 

If [lndepen[v, r] ==l&SCo[v, r] ==l&&Star[v, r] ==1, 
pset = KerAd[vi , b, n] ; 

If[r>l, Do [pset = Cap [pset, KerAd[v ± , : 

Length [pset] ^ 

Y= ^ f k ^psetlk, i]] b i; 

k=l 1 = 1 

While [Maximal [v, r, b, n] ==0, r = r + 1; Addone[r, v, b, i 
Print ["The rank is ", r, " and a Cartan subalgebr 

Do [comb = Solve [v q == X] ; 

Print [ (HoldForm[v] ) , " = ", 

Sum[ (a p // . comb[ [1] ] ) * (HoldForm[b] ) , {p, 
Print [False] ] ! 
By using the first element of the basis we have 



14 
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RankOf g2 [ {bi } ] 

The rank is 2, and a Cartan subalgebra: 

V] = bi 
v 2 = b u 

and the rank is 2. We can also try with 

RankOf g2 [{5 b 5 + bi } ] 

The rank is 2, and a Cartan subalgebra: 

V] = b\ + 5 bj 
v 2 = b\ + 5 b 9 



■ Lie algebra split g2 

We now deal with a particle physic model: the Lie algebra split g2 ■ A basis 
of this Lie algebra is 

BSplitg2 : = Module [ { } , 

n = 14; t = 8; e = Table[0, {i, 1, t}, { j, 1, t}] ; 
Do[m = 6; m|[p, gj = 1; F p , q = m, {q, 1, t}, {p, 1, t}] ; 
Do[J p , q = F p , q -F q , p ; Z p , q = F p , q +F q , p , {p, 1, t}, {q, 1, t}] 

bl = J2,3 - Jl,6 ; b2 = -J4,2 + J6,8 ', = %2,S + %4,1 } 

b4 = z 2 ,6 - z 4j8 ; bs = z 2 ,7 - z 4/5 ; b6 = z 2 ,s + Z4 ; 6 ; 

b7 = J3,4 + ^7,8 ; bs = Z 3/5 - Z 4/6 bg = Z 3;6 + Z 4;5 ; 
bio = Z 3/7 - Z 4/8 ; bu = Z 3j8 + Z 4;7 ; bl2 = - J6,5 _ J7,8 ; 
bl 3 = Js,7 + J6,8 ; bl4 = Js,8 + ^7,6 ', 

t t n 

Scala[x_, y_] :=^^x[[i, jj y[[i, jj ; X = ^] a± b ± ; ] 

i=l j=l i=l 

and we can make 



BSplitg2 

in order to have that a generic element is as follows: 
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'0 





























ai 


a 2 


a 3 


04 


as 


a 6 





—a\ 





a 7 


a 8 


ag 


aio 


an 





-a 2 


—a-j 





flg - as 


flg - <28 


a 3 + an 


-a 4 - «i 





a 3 


«8 


ag - as 





ai2 


a B 


fl!4 





CI4 


at) 


a 6 - a 8 


-«12 





a\ - ai4 


a 2 + an 





a 5 


flio 


a 3 + an 


-an 


^14 — a\ 





a 7 - an 


.0 


a 6 


a n 


-04 - aio 


-a u 


-a 2 -a [3 


a 12 - a 7 






We define then 



RankOf Splitg2 [dat_] : = Module [ { } , 

r = Length [dat] ; 

Do[v k =dat[[k]], {k, 1, r}]; 

If [lndepen[v, r] ==l&SCo[v, r] ==l&SStar[v, r] = = 

pset = KerAd[vi , b, n] ; 

If[r>l, Do [pset = Cap [pset, KerAd[vi, b, n] 
U, 2, r}]]; 

Length [pset] ^ 

Y = ^ fk Jpsetlk, i]] bi ; 

k=l 1 = 1 

While [Maximal [v, r, b, n] ==0, r = r + 1; Addone[r, v, b, i 
Print ["The rank is ", r, " and a Cartan subalgebr 

Do [comb = Solve [v q == X] ; 

Print [ (HoldForm[v] ) , " = ", 

Sum[(a p //. comb[[l]]) * (HoldForm [b] ) p , {p, 1, n}]], 
{q, 1, r}], Print [False] ] ] 

to check that the rank of this Lie algebra is also 2, as we have 



RankOf Splitg2 [ {bi + b 2 } ] 

The rank is 2, and a Cartan subalgebra: 

Vi = b\ + b 2 



v 2 =2b 4 -b l0 +3b u 



or 
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RankOf Splitg2 [ {bi } ] 

The rank is 2, and a Cartan subalgebra: 

V] = bi 

v 2 = b 4 +b l0 

Lie algebra so(4,4)©so(2,2) 

Finally, we present the Lie algebra (see the references) so(4,4)©so(2,2). 
This Lie algebra is a model we have created that explains all the quarks and 
their corresponding antiquarks. We define a basis by making 



e = Table[0, {i, 1, 12}, {j, 1, 12}]; 
Do[m = 6; m|[p, qj = 1; 

F p , q =m, {q, 1, 12}, {p, 1, 12}]; 

DO [ Jp,q = Fp ; q - Fq ;P ; Z p/ q = Fp ; q + Fq ; p , 

{p, 1, 12}, {q, 1, 12}]; 
bas={Zi,5, Zi /6 , Zi ;7 , Zi /8 , 

%2,5r Z 2 ,6, Z 2 ,7» Z 2 ,8, Z 3;S , 
Z 3 ,6r Z 3/7 , Z 38 , Z 4;5 , Z 4j6 , 

Z 4,7 f Z 4j8 , Ji ;2 , Jl,3 , 
Jl,4r i J2,4, J3,4j Js,6/ 

J5,7j J5,8f J6,l r ^6,8, 
Jl,8r J9,10» Jll,12r Z 9; h, 

Z9,12, Zio,n, Zio,12}; 

n = Length [bas] ; 

t = 12; Do[b ± =bas[[i]], {i, 1, n}]; 

t t 

Scala[x_, y_] :=^^x[[i, jj y|[i, jj , 

i=i j=i 

n 

x = b±; 

i=l 

and then a generic element is 
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We define then 

RankOfModel [dat_] : = Module [ { } , 

r = Length [dat] ; 

Do[v k =dat[[k]], {k, 1, r}]; 

If [lndepen[v, r] ==1&&Co[v, r] ==l&SStar[v, r] = = 
pset = KerAd[vi , b, n] ; 

If[r>l, Do [pset = Cap [pset, KerAd[v ± , : 

Length [pset] ^ 

Y = ^ f k ^ pset [Ik, ij bi ; 

While [Maximal [v, r, b, n] ==0, r = r + 1; Addone[r, v, b, 1 
Print ["The rank is ", r, " and a Cartan subalgebr 

Do [comb = Solve [v q == X] ; 

n 

Print["v" q , " = ", ^ ( a P //• combEH) *b p ], {q, 1, 1 

Print [False] ] ] 

We can now compute the rank of this Lie algebra by using two elements: 



Pablo Alberca -Cdndido Martin 



RankOfModel [ {bi , b 6 }] 

The rank is 6, and a Cartan subalgebra: 

vi =bi 

v 2 = b 6 

v 3 = b 34 

v 4 = b 3l 

v 5 = b l6 

v 6 = fen 

or with three elements 

RankOfModel [ {bi , b 6 , bn } ] 

The rank is 6, and a Cartan subalgebra: 
V] =b x 

v 2 = b 6 
v 3 = 6 U 
v 4 = b 34 
v 5 = 6 31 
v 6 = 616 

This is the reason why this model explains six observables by acting on R 
(6 quarks + 6 antiquarks). 



References 

Alberca, P. and Martin Gonzalez, C. Automorphisms groups of composition algebras and 
quarks models. Hadronic Journal. Florida. 1997. 



Rank of a matrix Lie algebra 



Alberca, P. and Martin Gonzalez, C. Cdlculo efectivo del rango de las algebras de Lie de 
matrices. EACA. 1997. Universidad de Granada. 
Jacobson, N. Lie algebras. Dover Publications, Inc. 1979. 

Kaufman, S. Mathematica as a tool. An introduction with practical examples. Birkauser. 
1994. 

Postnikov, M. Lecons de geometrie: groupes et algebres de Lie. Mir. 1990. 

Wolfram, S. Mathematica. A system for doing mathematics by computer. Addison-Wesley. 
1988. 



